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We report a comprehensive experimental study and detailed model analysis of the terahertz dielec- 
tric response and density kinetics of excitons and unbound electron-hole pairs in GaAs quantum 
wells. A compact expression is given, in absolute units, for the complex-valued terahertz dielectric 
function of intra- excitonic transitions between the Is and higher-energy exciton and continuum lev- 
els. It closely describes the terahertz spectra of resonantly generated excitons. Exciton ionization 
and formation are further explored, where the terahertz response exhibits both intra-excitonic and 
Drude features. Utilizing a two-component dielectric function, we derive the underlying exciton 
and unbound pair densities. In the ionized state, excellent agreement is found with the Saha ther- 
modynamic equilibrium, which provides experimental verification of the two-component analysis 
and density scaling. During exciton formation, in turn, the pair kinetics is quantitatively described 
by a Saha equilibrium that follows the carrier cooling dynamics. The terahertz-derived kinetics 
is, moreover, consistent with time-resolved luminescence measured for comparison. Our study es- 
tablishes a basis for tracking pair densities via transient terahertz spectroscopy of photoexcited 
quasi-two-dimensional electron-hole gases. 

PACS numbers: 78.47.-p, 73.20.Mf, 78.67.De 



I. INTRODUCTION 

The terahertz (THz) frequency electromagnetic response 
provides important insight into low-energy excitations 
and many-body correlations in condensed matter^ In 
semiconductors, Coulomb interactions lead to the forma- 
tion of excitons from unbound electron-hole (e-h) pairs. 
Microscopic interactions on ultrashort time scales deter- 
mine the dynamics of energy relaxation, dephasing, diffu- 
sion, or species inter-conversion of excitons and unbound 
pairs. Besides opto-electronic applications, understand- 
ing the low-energy structure and dynamics of these quasi- 
particles is fundamentally important, e.g., for the explo- 
ration of low-temperature collective phenomena^i 3 -^ 

Excitons were extensively studied via optical absorp- 
tion, photoluminescence (PL), or nonlinear experiments 
at the semiconductor band gap. 6 Despite their success, 
these techniques rely on e-h pair creation and annihila- 
tion which entails important limitations. In particular, 
momentum conservation often restricts the sensitivity of 
PL to a subset of excitons around center-of-mass momen- 
tum K ~ 0. The intensity then depends not only on the 
density of the exciton gas but also on its detailed distribu- 
tion function. Higher- .fT excitons can be detected in ma- 
terials with strong electron-phonon coupling, giving in- 
sight, e.g., into exciton thermalization and formationj^^ 
However, the determination of absolute densities is ex- 
ceedingly difficult lacking precise knowledge of the col- 
lection efficiency and inter-band dipole moment. 

In contrast, intra-excitonic THz transitions between 
the low-energy internal levels of excitons represent a fun- 
damentally different tool. They measure the coupling 



between the Is exciton ground state and higher relative- 
momentum states, detecting excitons largely indepen- 
dent of K . Intra-excitonic probes can determine abso- 
lute exciton densities, since they measure - analogous to 
atomic absorption spectroscopy - existing excitons with 
predictable cross section. These THz transitions are inde- 
pendent of the interband dipole moment, which renders 
them suitable probes in the search for exciton condensa- 
tion. Importantly, THz fields are also equally sensitive 
to unbound e-h pairs (free carriers). 

While initial studies were scarce, intra-excitonic spec- 
troscopy recently emerged as a powerful tool to in- 
vestigate the low-energy resonances and dynamics of 

c::c jt on:: 10.11.12. 13.14.15.16, 17.18.19.20, 21.22.23.24.25.26.27.28 f u _ 

eled in part by rapid advances in ultrafast THz 
technology. 2 ^ Optical-pump THz-probe experiments are 
of particular interest, as they yield both real and imagi- 
nary parts of the transient response functions. Notably, 
we reported distinct THz signatures of conducting and in- 
sulating phases during exciton formation and ionization 
in GaAs quantum wells (QWs)fi 7 - Such complex- valued 
THz spectra place strict boundaries on theoretical mod- 
els, providing exceptional potential for further analysis 
to directly determine exciton and free carrier densities. 

Here, we report a detailed model analysis and exper- 
imental study of the transient THz spectra and density 
kinetics of quasi-two-dimensional (2D) excitons and un- 
bound e-h pairs in GaAs QWs. An accurate evaluation 
procedure for the optical-pump THz-probe signals is de- 
scribed for the multi-layer QW geometry. A key aspect 
of this work is the derivation of the intra-excitonic dielec- 
tric function scaled in absolute units, as a compact ex- 
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pression to determine pair densities from measured THz 
spectra. It closely describes both shape and amplitude of 
the THz response during population decay of resonantly- 
generated excitons. In contrast, the THz line shapes dur- 
ing exciton ionization and formation exhibit both intra- 
excitonic and Drude-like features. A two-component di- 
electric function is applied to infer the underlying exciton 
and free-carrier fractions. After ionization, this reveals 
a mixture that agrees closely with the thermodynamic 
equilibrium (Saha) predictions - thus experimentally ver- 
ifying the two-component analysis. Its application to ex- 
citon formation reveals a surprisingly simple yet quanti- 
tative description of the dynamics, when cooling of the e- 
h gas is taken into account. This THz-derived scenario is 
found to consistently describe the luminescence dynamics 
measured for comparison. Hence, this paper establishes 
a basis for gauging exciton and free-carrier densities in 
the THz dielectric response of quasi-2D e-h gases, and 
demonstrates its application to tracing the time evolu- 
tion of multi-component phases during exciton ionization 
and formation. 

In the following, Sec. |TT] explains the experimental 
methods, while the THz response of resonantly-generated 
excitons and unbound pairs is discussed in Sec. IIIII Mod- 
els of the intra-excitonic and Drude response are pre- 
sented that enable determination of absolute pair densi- 
ties. Section [IV] discusses the THz response and kinet- 
ics during exciton ionization at elevated lattice tempera- 
tures. In Sec. El exciton formation after non-resonant 
excitation is analyzed and compared to PL kinetics. 
Appendixes FA] and [B] derive the complex transmission 
function of the quantum-well structure, and the intra- 
excitonic dielectric function, respectively. 



II. EXPERIMENTAL TECHNIQUE 

In the experiments, we utilize optical-pump THz-probe 
spectroscopy to investigate transient changes in the THz 
conductivity of photoexcited e-h gases in GaAs QWs. 
The high-quality sample studied here consists of a stack 
of ten 14-nm wide, undoped GaAs wells separated by 
10-nm wide Alo.3Gao.7As barriers, grown via molecu- 
lar beam epitaxy on a GaAs substrate^ The structure 
is embedded in 500 nm-thick Alo.3Gao.7As spacer lay- 
ers. To avoid THz absorption of photoexcited carriers 
in the substrate, it was removed by selective etching 
after attaching the QW side to a 0.5-mm thick (100) 
MgO substrate.™ The corresponding low-temperature 
near-IR absorption spectrum is shown by the solid line 
in Fig. [2(a). The Is heavy-hole (HH) exciton line at 
1.540 eV dominates the spectrum, with a line width of 
0.8 mcV [full width at half maximum (FWHM)]. The 
2s HH and Is light-hole (LH) absorption lines, and the 
inter-band continuum follow at higher photon energies. 

The experimental setup is shown in Fig. [2(b). We 
utilize amplified pulses at high repetition-rate to allow 
for both sensitive THz detection and sufficiently intense 
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FIG. 1: (color online), (a) Near-IR absorption of the GaAs 
multi-quantum- well sample at T = 6 K (solid line) . A typical 
laser spectrum is shown unshaped (dotted line), and shaped 
for resonant ls-HH excitation (dashed line) or continuum ex- 
citation (dashed-dotted line), (b) Experimental setup for 
optical-pump THz-probe spectroscopy. The THz section is 
purged with dry nitrogen to avoid far-infrared absorption in 
air. 



pulses to excite the intrinsically large THz probe area. At 
the outset, a 250-kHz Tksapphirc amplifier system (Co- 
herent RegA) delivers 150-fs near-IR pulses at 800 nm 
wavelength. A fraction of the output is used to gen- 
erate THz probe pulses via optical rectification and to 
detect them via electro-optic sampling, each in a 500-^m 
thick (110) ZnTe crystal^ The THz pulses span the 2- 
12 meV (« 0.5-3 THz) spectral range. The THz beam 
is recollimated and focused with off-axis parabolic mir- 
rors onto the sample mounted in a cold-finger cryostat. 
For optimal time resolution and to avoid spectral pump- 
probe artifacts, the THz field is scanned via the genera- 
tion pulse delay and the pulses are incident from the sub- 
strate side^ The THz focus size is frequency-dependent, 
with a FWHM diameter of around 1 mm at 1.5 THz.— 
We employ a 2-mm diameter aperture to limit the THz 
probe to the photoexcited region. 

A second part of the near-IR amplifier output is used 
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FIG. 2: (color online), [(a)-(c)] THz response at XL = 300 K, at delay At = 10 ps after non-resonant excitation 100 meV above 
the band gap with excitation fluence F = 0.6 ^tj/cm 2 . The THz transients in (a) are the reference (dashed line) and pump- 
induced change (solid line, multiplied x50). Panels (b) and (c) show corresponding spectra Atri and Aei from experiment 
(dots), and a Drude model (solid line) with n e h — 2 x 10 10 cm -2 and Fd = 4.8 meV. [(d)-(f)] THz response for resonant 
ls-HH excitation at TL = 6 K, with F = 0.14 /xj/cm 2 and At = 10 ps. Solid lines in panels (e) and (f): intra-excitonic model, 
tlx = 2.7 x 10 10 cm~ 2 . (g) Intra-excitonic model as in panel (e), but Fhb = Tb c = 1 /UeV. Solid line: ls-np transitions; dashed 
line: ls-continuum. 



for photoexcitation. The full laser spectrum is shown as 
the dotted line in Fig. HJa). For selective excitation of 
excitons or unbound pairs, the pulses must be spectrally 
shaped. This is achieved in a reflective, zero-dispersion 
stretcher consisting of a 1200 1/mm grating, 200- mm fo- 
cal length lens, and an adjustable slit. Pump spectra are 
narrowed to w 1-2 meV width, with typical curves shown 
as the dashed and dashed-dotted curves in Fig.QJa). This 
yields an overall 1-2 ps time resolution. On-line spectral 
characterization of the pump light transmitted through 
the sample allows for precise determination of the exciton 
line position, for optimal spectral overlap and resonant 
excitation at varying sample temperatures. 

The THz studies are complemented by time-resolved 
PL on the identical sample. Such data were taken with 
the sample mounted in a vapor-flow cryostat and pho- 
toexcited by a 76-MHz Tusapphire oscillator, spectrally 
narrowed using a 1-meV wide interference filter. Linearly 
polarized excitation was used to minimize spin effects. 
The luminescence was temporally resolved with a Hama- 
matsu streak camera. These measurements necessitated 
a smaller pump spot diameter (140/^m), which however 
remains well beyond estimated carrier diffusion lengths 
within the 1-ns time window. 

In the THz experiments, the dynamics of the dielectric 
response is determined as follows. After photoexcitation, 
the THz dielectric function of the QW layers in equilib- 
rium, denoted by e(w), transiently changes to the modi- 
fied value e(w) + Ae(w). Thus, the induced change Ae(u) 
must be determined for each fixed time delay At between 
the arrival of pump and probe pulses on the sample. 
For this, we measure the THz reference probe field E(t) 



(transmitted through the sample in equilibrium) and its 
pump- induced change AE(t). Typical time-domain THz 
signals are shown in Figs. [2la) and[2]^d). Fourier trans- 
formation provides the corresponding frequency-domain 
fields E(ui) and AE(u). The change in the dielectric re- 
sponse is then obtained as a function of these fields, 

which depends on the sample geometry. An analyti- 
cal expression / that takes into account the multilayer 
structure of our QW sample is derived in Appendix [X] 
[Eq. (lATOl) ]. 

The current response J(w) = a{uj)E(uj) of the many- 
particle system to the incident transverse electromag- 
netic field is given by the optical conductivity (t(oj) = 
ai(u) + ia^ioS)- It is connected to the dielectric func- 
tion via u{lj) = iujeo[l — e(u>)]. In the following, we will 
express the transient THz response as 

Ae(w) = Aei(w) + — A(7i(w). (2) 
e uj 

Here, the induced conductivity Aa±(uj) is a measure of 
the absorbed power density and allows for analysis of os- 
cillator strengths. The dielectric function change Aei(w), 
in turn, provides a measure of the inductive, out-of-phasc 
response. As evident below, the availability of both Ao^ 
and Aei is a key to distinguishing different contributions 
to the multi-component THz spectra. 
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III. TERAHERTZ RESPONSE OF E-H PAIRS 

We first discuss experiments that probe transient changes 
in the THz dielectric response of unbound e-h pairs after 
non-resonant photoexcitation into the band-to-band con- 
tinuum. Data are shown in Figs.[2ja)-[2](c) for lattice tem- 
perature Xl — 300 K, which ensures rapid ionization of 
the e-h pairs. As evident in the time-domain THz traces 
in Fig. HJa), the pump-induced field change (solid line) 
resembles the reference (dashed line) with a phase shift, 
pointing to a spectrally broadband response. This is con- 
firmed by AfTi(ci;) and Aei(u>) shown as dots in Figs. [21(b) 
and^c). The large low- frequency conductivity Aeri un- 
derscores the conducting nature of the unbound pairs, 
while the dispersive Ae x < is characteristic of a zero- 
frequency, Drude-like oscillator. Indeed, the response is 
well described by the Drude model (solid lines) 



Ae(w) = n eh Ae D (w) = n ch 



dweofi (w 2 + iLoTv) 



, (3) 



where n e h is the e-h pair sheet density per well, c?w is 
the QW width, Td is the Drude scattering rate, and 
fx = ra e mh/(m e + Wh) is the reduced mass (m c , mi, 
are e and h effective masses, respectively). The den- 
sity n e h from the Drude fit closely agrees with the value 
(1.9 x 10 10 cm~ 2 ) estimated from the excitation fluence 
and sample parameters. 

Next, we will discuss the transient THz response af- 
ter resonant excitation at the ls-HH exciton line, for 
7l = 6 K. The measured THz field change in Fig.OJd) ex- 
hibits a complex shape with time-dependent phase shift. 
Figures [He) and[^f) show the corresponding THz spec- 
tra (dots). Here, the conductivity Acti is characterized 



by a distinct, asymmetric peak around Hlo 7 meV. The 
dielectric function change Aei, in turn, shows an oscil- 
latory response around the same photon energy. This 
represents a new low-energy oscillator, absent in equilib- 
rium, which can be explained by transitions between an 
exciton's internal degrees of freedom. The 7 meV con- 
ductivity peak arises from the Is — > 2p transition between 
the exciton levels, in concordance with GaAs/AlGaAs 
QW exciton binding energies 3 -^ and the Is — 2s splitting 
in Fig. [IJa). The vanishing low- frequency conductivity 
Aei of the intra-excitonic response is a signature of the 
insulating nature of the charge- neutral excitons. Com- 
pared to the Drude response, the opposite sign of Aei at 
low frequencies enables further discrimination between 
the THz response of excitons and unbound pairs. 

For a quantitative description we have performed cal- 
culations, scaled in absolute units, of the intra-excitonic 
contribution to the dielectric function. The model takes 
into account 2D bound and continuum hydrogenic wave 
functions, where the Bohr radius 



and binding energies 



A-nh e e s 
e 2 /xA 



E n — 



8"7reoe s a 



(n+1/2)" 



(4) 



(5) 



are scaled by the reduced mass /j, and static dielectric con- 
stant e s , and by a parameter A which scales the Coulomb 
potential to take into account the finite well width. 36 
The complex- valued dielectric THz response from intra- 
excitonic transitions between Is and higher bound and 
continuum states is then given by 



o„2„2 °° Tp _ Tp i 1 

AeM = n x A ex ( W ) = nxf^-{£ ( . E g . - ■ (1 - -^) 2 " ^ 3 + (6) 

dwK eo ~[ ([ "ft ° ] 2 ) - wF bb 2n + 2 n + 1 

E + E(k) o?k (^)-a dk} 



{[ E 2± E^ ]2 _ UJ 2 ) _ lujThc i + e-S (l + {ak/2}*y 



where fix is the Is exciton sheet density, Tbb and 
Tbc are level broadenings for bound-bound and bound- 
continuum transitions, and E(k) = h 2 k 2 /2n is the kinetic 
energy of the continuum states. Details of the derivation 
are given in Appendix [B] To illustrate the underlying 
transitions and oscillator strengths, Fig. [2^g) shows a 
calculated intra-excitonic line shape for (unrealistically 
low) 1 /ieV broadening. It consists of the ls-2p peak 
and transitions into higher bound np levels (solid line), 
and transitions into the continuum of unbound e-h pairs 
(dashed). As for three-dimensional intra-atomic tran- 
sitions, bound-bound and bound-continuum transitions 



match up smoothly at the continuum edge. 

We now compare this model with the experimental 
data. The calculated intra-excitonic response is shown 
as solid lines in Fig. ^e) and H^f), which reproduces 
the shape of the experimental data extremely well. The 
sharp lines are now absent due to realistic broadening, 
and the response is dominated by the ls-2p peak and 
higher-energy shoulder. Most of the model parameters 
are severely restricted: for the GaAs QWs e s = 13.2 
and /j, = 0.054m ( TO o : f ree electron mass), using m c = 
0.0665m and — 0.28to averaged in /c-space over 
w 1/a to account for the exciton wave function^ More- 
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FIG. 3: (color online), (a) Transient THz response (dots) 
for different delays At, after resonant ls-HH excitation as 
in Figs. [2|e) and[2|f). Lines: intra-excitonic dielectric model, 
with r bb = 0.9 meV (0 ps), 0.8 meV (5 ps), 1.03 meV (30 ps), 
1.1 meV (300 ps), and 0.9 meV (1000 ps). Curves shifted 
vertically are scaled identically, (b) Exciton density nx (dia- 
monds) from the model in panel (a). Dashed line: exponential 
decay with r = 1087 ps. Solid line: PL intensity after reso- 
nant excitation, for n ~ 2 x 10 10 cm" 2 and TL = 10 K. 



over, A = 0.678 must be chosen to reproduce the observed 
ls-2p level spacing. This leaves, as the only free param- 
eters, the Is exciton density nx and level broadenings 
r. The best agreement in shape in Figs. [2^e) and [2jf ) 
is obtained with Tbb = 0.8 meV for bound-bound tran- 
sitions and Tbc = 2.2 meV for bound-continuum tran- 
sitions. The larger Tb c reflects increased scattering of 
continuum final states, as corroborated by a compara- 
ble Drude width of the non-resonantly excited, T = 6 K 
response discussed further below. For the experimental 
data in this paper, the consistently best fit was obtained 
by keeping Tb c fixed while varying the intra-excitonic 
broadening Tbb- The Is exciton density from the model 
fit (n x = 2.7 x 10 10 cm- 2 ) in Figs. H» and^f) com- 
pares well with the density 2.1 x 10 10 cm" 2 estimated 
from the pump flux (0.14 /LtJ/cm 2 ) in the experiment af- 
ter accounting for sample absorption, spectral overlap, 
and reflection losses of the cryostat windows and sample. 
Work at densities well below the present 10 10 cm -2 range 
is desirable, but will necessitate further improvements in 
the sensitivity of the THz measurement technique. 

With respect to absolute density scaling, it should 
be emphasized that both the Drude and intra-excitonic 
models above fulfill the "partial oscillator strength sum 



rule" 

/ <Tl{U3)du=-- — , (7) 

as expected in a parabolic band approximation. The to- 
tal photoexcited sheet density n of bound and unbound 
pairs can thus be directly obtained from the integral of 
A(Ti (to below the onset of inter-band transitions), which 
corresponds to the induced intra-band spectral weight. 
This underscores the capability of THz spectroscopy to 
determine, unlike luminescence, absolute densities of ex- 
citons and unbound e-h pairs both at K = and outside 
the optically-accessible momentum range. 

Transient THz spectra at several different pump-probe 
delays are shown in Fig.[3ja). Initially, a coherent Is ex- 
citon polarization is created by the near-IR pump pulse 
which dephases within a few picoseconds into an incoher- 
ent exciton population^ In the coherent regime directly 
after excitation, the THz response deviates noticeably 
from the model line shape, as evident at At = ps in 
Fig. EJa) . At later times (At > 5 ps) the intra-excitonic 
dielectric function well describes the transient THz re- 
sponse. The conductivity decays in amplitude but retains 
its peaked line shape, evidencing directly the decay of ex- 
citonic populations. Changes in the broadening Tbb are 
minor, and may result from a time-varying temperature 
of the exciton gas due to recombination heating. Densi- 
ties derived from the model fits are charted in Fig. [3jb) 
(diamonds). An exciton recombination time r ~ 1 ns 
is obtained from the single-exponential fit (dashed line). 
The radiative decay of free K = excitons in QWs is pre- 
dicted to be as short as 10 ps in the idealized case with- 
out any dephasing. PL experiments on a highest-quality 
single-QW indeed exhibited decays as short as 40 ps, for 
narrow well width (4.5 nm), low density (3 x 10 9 cm" 2 ) 
and very low temperature (1.7 K)i^s However, for com- 
parison with the THz-derived dynamics we must consider 
the case of wide QWs and densities exceeding 10 10 cm" 2 , 
where PL decay times are much longer and consistent 
with our result i 39 ' 40 

We also compare the THz results to our time-resolved 
PL measurements after resonant ls-HH excitation, in- 
dicated by the solid line in Fig. \3\b). The PL ex- 
hibits an initially faster decay, which can arise from un- 
suppressed spin relaxation effects, coherent emission, or 
pump scattering, effects that typically complicate lu- 
minescence transients. At later times the PL agrees 
well with the THz dynamics. Unlike the THz exper- 
iment, however, absolute densities are intrinsically dif- 
ficult to obtain from the PL intensity, despite that in 
this resonantly-excited case it reflects recombination of 
mainly cold excitons in K « luminescent states. 

IV. EXCITON IONIZATION 

At increased lattice temperature, the temporal dynamics 
and shape of the transient THz spectra undergo exten- 
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FIG. 4: (color online). Transient THz spectra at lattice temperatures (a) Tl = 20 K, (b) 40 K, and (c) 80 K after resonant ls- 
HH excitation with fluence of 0.14 /iJ/cm 2 . Dots: induced conductivity and dielectric function change at indicated pump-probe 
delays At. Curves are shifted vertically, but scaled equally. Solid lines: two-component model; dashed lines: intra-excitonic 
model only, with time-varying density and broadening^! 



sive changes. Figure [4] shows the response after resonant 
HH excitation, for three representative temperatures. At 
20 K, Af7i(w) is well described by a sharp exciton line 
shape directly after excitation, but broadens noticeably 
with increasing time delay [Fig. H^a)]. Simultaneously, 
Aei(w) flattens out. As shown in Figs. (Hb) andQJc), 
these changes occur faster and become even more en- 
hanced as the lattice temperature is further elevated to 
40 and 80 K. Two important hallmarks of a Drude re- 
sponse appear: at the low-frequency end of the spectrum 
significant conductivity Aei builds up, and the induced 
dielectric function Aei increasingly deviates towards neg- 
ative values. The dynamics in Fig. 3] thus evidences the 
generation of unbound e-h pairs, which indicates thermal 
ionization of the resonantly excited HH excitons. 

A. Two-component analysis 

To obtain a quantitative picture of the ionization pro- 
cess, we need to describe the complex spectra at all delay 
times. The intra-excitonic model function alone (dashed 
lines in Fig. 0| is clearly insufficient. Indeed, the intra- 
excitonic Aei always remains positive below the ls-2p 
oscillator frequency (y < 1.7 THz), regardless of the 
amount of broadening. This underscores the importance 
of measuring both real and imaginary parts of the THz 
response. In order to take into account the simultaneous 
existence of excitons and unbound e-h pairs, we imple- 
ment a two-component dielectric function 

Ae(w) = n x Ae x (w) + n hAe D (w) , (8) 

where Aex is the intra-excitonic dielectric function and 
Aed is the Drude response described above. This model 
was fitted to each THz spectrum by varying the densities 



^ehj n x and the broadening parameters. Resulting model 
functions are shown as solid lines in Fig. 2J They describe 
the experimental data well, while neither an exciton nor 
a Drude model alone can account for the response at long 
delay times and elevated temperatures. The fit param- 
eters are strongly constrained due to important spectral 
differences between the response of excitons and unbound 
e-h pairs, and by the need to explain both Aai(u) and 
Aei (a;) simultaneously and over a broad spectral range. 
These THz spectra are reproduced with a fixed binding 
energy in the model; we verified that, using the equa- 
tions in Ref. l42l . that renormalization due to free-carrier 
screening should indeed remain < 10% for our condi- 
tions. However, a distinct time-dependent broadening is 
observed which reflects changes in homogeneous dephas- 
ing of both the Is and 2p levels due to exciton-free car- 
rier and exciton-phonon interactions, with possibly inho- 
mogeneous contributions due to a momentum-dependent 
binding energy^ 

With knowledge of n cn and nx in absolute units, we 
can obtain the exciton fraction 

/x^^^, (9) 

as a measure of the admixture of bound e-h pairs to 
the many-particle system. Figure E]Ja) shows the tem- 
poral dynamics of fx obtained from two-component fits 
to the experimental data. At higher lattice temper- 
atures, fx decays with time until it reaches a quasi- 
equilibrium value. As evident, with rising lattice temper- 
ature the ionization becomes faster while the residual, 
quasi-equilibrium exciton fraction at long delay times 
decreases. The total pair density is fairly constant for 
all temperatures, as shown in Fig. EJb), and decays 
monotonously with time. 

Figure [U^c) shows the temperature dependent ioniza- 
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FIG. 5: (color online), (a) Ionization kinetics after resonant 
excitation at indicated temperatures. Symbols: exciton frac- 
tion fx from the two-component analysis. Lines: fit with 
exponential decay plus offset, (b) Total pair density N for At 
= 5 ps (dots), 40 ps (squares), and 300 ps (diamonds), (c) Ex- 
citon ionization rate (dots) from the initial decay of fx com- 
pared to LO-phonon scattering (line) with Flo = 30 THz. (d) 
Exciton fraction fx at long times (squares), compared to the 
Saha model for N = 2 x 10 10 cm -2 (solid line), 1 x 10 10 cm -2 
(dashed line), and 4 x 10 10 cm -2 (short-dashed line). 

tion rate, as derived from the initial decay rate of the 
exciton fraction, dfx/dt\ t =o from Fig. Efa). We can ex- 
plain the large temperature dependence above XL > 50 K 
by ionization through LO phonon absorption. The scat- 
tering rate of this process is given by T = IYo ■ n{Ti), 
where n(T) = (exp(£lLo / k^T) — is the Bose oc- 
cupation with SIlo = 36.6 meV for GaAs. To de- 
scribe the experiment, a corresponding model function 
with Flo = 30 THz is shown in Fig. |Sfc) (solid line). 
This value for Tlo is in excellent agreement with the 
LO phonon scattering rates of excitons derived from the 
broadening of the near-IR exciton absorption finest 



where N = n h + «x is the total e-h pair density. The 
reduced mass \i = 0.054m and binding energy E Q — 
7.7 meV are retained from the above lineshape model. 
For a given total pair density N, Eq. [TU] then yields the 
temperature dependence of nx, and hence fx{T), in the 
thermodynamic equilibrium. 

Figure (Hd) compares the Saha prediction of fx(T) 
(lines) with the experimentally-derived values (squares) 
for long delay times. A pair density N = 2 x 10 10 cm -2 
(solid line) yields surprisingly close agreement with 
the temperature dependence of the experimentally- 
determined exciton fraction. This extends previous PL- 
based studies of the 100-300 K range (Ref. H^) to tem- 
peratures below 100 K and provides quantitative density 
information accessible only to THz probes. Importantly, 
not only the shape fx(T) but also the absolute density 
N underlying this Saha model curve agrees well with the 
values obtained from the THz spectra, as evident from 
the total pair density at long delay times (At = 300 ps) 
shown as diamonds in Fig. [5jb). 

An alternate calculation of the THz response of mixed 
e-h gases was also reported^ which for THz spectra with 
overwhelming 90% intra-excitonic oscillator strength pre- 
dicts an exciton fraction of only 10%. 48 Consider this pic- 
ture applied, e.g., to our THz spectra after ionization at 
T = 20 K corresponding to fx = 85% [Figs. HJa) and 
\§[d)] . The model of Ref. |48| accordingly implies an exci- 
ton fraction < 10%, more than six times the ionization 
level possible by thermal excitations at the given tem- 
perature, density, and binding energy. Also, since ls-2p 
THz absorption directly measures exciton populations, 16 
the exciton density is ~ 2 x 10 10 cm~ 2 [Fig.^a), 400 ps] 
using our gauge from Sec. Mil A 10% exciton fraction 
would then imply a total pair density of 2 x 10 11 cm -2 , 
exceeding the known, absorbed photon density more than 
seven times. Given these stark discrepancies, we cannot 
employ the model of Refs. [47] and|H for quantitative in- 
sight into THz spectra of mixed e-h gases. 

In contrast, the self-consistent, quantitative agreement 
shown in Fig. [5] between (i) the Saha-model densities 
and exciton fraction founded on basic thermodynamic 
relations, and (ii) the experimentally-derived total pair 
density and exciton fraction from our analysis during the 
quasi-equilibrium at long delay times provides a clear and 
direct validation of our two-component dielectric function 
analysis of the transient THz spectra. 



B. Thermodynamic quasi-equilibrium 

We can now compare the quasi-equilibrium observed 
at long delay times to the predictions of the so-called 
Saha equation. The latter describes the densities of ex- 
citons and free carriers after statistical equilibration of 
their chemical potentials in the Boltzmann limit. For a 
2D gas of e-h pairs, it reads i 44 i 45 ' 46 

{N - nx? = M e-^, (10) 

nx 2nh 



V. EXCITON FORMATION 

Having verified the applicability and correct gauge of 
the above dielectric function for the study of mixed e-h 
gases, we can now analyze the kinetics of exciton forma- 
tion from unbound pairs. For this, the QWs are excited 
non-resonantly above the band gap. Transient THz spec- 
tra are shown in Fig. [6l They exhibit a complex dynam- 
ics which evolves from a broad Drude-like response into 
an intra-excitonic line shape. Directly after excitation, 
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FIG. 6: (color online). Transient THz spectra (dots) at 
Tl = 6 K after non-resonant excitation into the continuum 
(at 1.561 eV with fiuence 0.2 /xJ/cm 2 ). Solid lines: two- 
component model; dashed lines: intra-excitonic model only, 
with broadening Vbb = 3 meV (5 ps), 2.8 meV (60 ps), 2 meV 
(300 ps), 1.5 meV (600 ps), and 1 meV (1100 ps) and densities 
as in Figs.EIa) and^b). 

the low- frequency conductivity in Aeri and the —1/uj 2 
dispersion in Aei reveal a conducting, Drude-like phase. 
However, a broad excitonic peak is also evident in Aa\ at 
these early delay times, rendering the spectra similar to 
the mixed ionized phase discussed above. With increas- 
ing delay the Drude component decays, while spectral 
weight builds up around the \s-2p transition until the 
intra-excitonic line shape is restored. Exciton formation 
thus proceeds on two different time scales: fast formation 
of a large exciton fraction directly after excitation, and 
a slow binding of remaining free carriers into excitons 
within several 100 ps. 

A. Formation kinetics 

For quantitative insight, we analyzed the transient 
spectra with the two-component dielectric function (solid 
lines in Fig. [6J. The corresponding exciton and un- 
bound pair densities and the exciton fraction are shown 
in Figs.[7]Ja)-|7fc), as obtained from the above data (solid 
dots) and from a data set with ~ 3.5 times higher den- 
sity (open circles). In both cases, the analysis confirms 
a fast initial generation of an appreciable exciton frac- 
tion fx ~ 40%, followed by slower transfer of unbound 
e-h pairs into excitons and an eventual decay of the exci- 
ton density due to recombination. At the longest delays 
(At w 1 ns) almost all carriers are bound into excitons 
with fx « 90%. 

After non-resonant excitation, photoexcited e-h gases 
are known to thermalize on a 100-fs time scale into a 
Fermi distribution, which cools to the lattice tempera- 
ture via emission of optical and acoustic phonons. 6 This 



raises the question whether Figs. [7ja)-[7fc) reflect a pair 
kinetics where (i) exciton formation is much faster than 
cooling, such that the carrier gas maintains thermody- 
namic equilibrium at each time, or where (ii) the time 
evolution is limited by slower, bi-molecular e-h pairing in- 
teractions resulting in a measurable, non-equilibrium de- 
viation from the Saha condition. To resolve this matter, 
we have calculated the cooling dynamics of the quasi- 
2D carrier temperature Tc(t) by integrating the time- 
dependent change 

dT c = (dE/dt) Tc (dE/dt)^ 
dt 2k B a 2k B a ' ^ ' 

where the second term ensures equilibration at the lat- 
tice temperature. Here, (dE/dt)r is the energy-loss rate 
of the e-h gas at temperature T, for which we directly 
employed the values obtained for GaAs QWs by Leo 
et alA^ Moreover, a is the well-known reduction fac- 
tor due to hot-phonon effects, and a = 3 was chosen 
in agreement with previous worki 49 i 50 i 51 Given the exci- 
tation conditions, we estimate an initial carrier temper- 
ature T c (0) = 69 K after distributing the excess energy 
of the photoexcited carriers equally among electrons and 
holes^i The resulting cooling curve is shown in Fig.[7]^d). 
It is characterized by a quick drop to a temperature of 
~ 40 K due to emission of LO phonons, after which the 
cooling proceeds via acoustic phonons on a much longer 
time scale. The values and overall shape of this cooling 
dynamics closely agrees with previous work i 46 i 50 

We can now compare the pair densities and exciton 
fraction obtained from the transient THz spectra with 
the predictions of the Saha equation. This model de- 
scription assumes that the excitons and free carriers 
form a quasi-equilibrium that conforms to the time- 
dependent carrier temperature Tc(t), via equilibration 
processes that are fast compared to the overall dynam- 
ics. Accordingly, the Saha equilibrium exciton frac- 
tion and densities were calculated from Eq. [TUl with 
N = nx(t) + fieh(i) a s obtained from the measured THz 
spectra and with T = Tc(t). The resulting curves are 
shown in Figs.[7Ja)-[7|c) as solid and dashed lines, clearly 
yielding an extremely close quantitative description of 
the experimentally-derived densities and exciton fraction. 
Hence, the thermodynamic equilibrium model provides 
for a surprisingly simple description of the formation ki- 
netics. 

The capability of the Saha quasi-equilibrium to explain 
the time-dependent exciton formation kinetics agrees 
with observations in several PL studies j 46 i 52 However, 
it should be emphasized that our observations do not 
rule out the influence of bi-molecular interactions in the 
formation kineticsj 50 ' 53 ! 54 Indeed, as explained by De- 
veaud et al.,^ in the 10 10 cm~ 2 density range appli- 
cable to our current study the Saha and bi-molecular 
rate-equation models are expected to be largely com- 
mensurate, while they deviate at lower densities (see 
Fig. 7 of Ref . l55t) . This motivates future work with even 
more sensitive THz probes to explore formation in the 
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FIG. 7: (color online). Pair density and luminescence during 
exciton formation, [(a)-(c)] Exciton density, free-carrier den- 
sity, and exciton fraction from the two-component model in 
Fig. [5] (solid dots) and from data at ~ 3.5 x higher density 
(open circles). The Saha model (see text) is shown as solid 
and dashed lines for low and high densities, respectively, (d) 
Simulated carrier temperature Tc for a = 3, TL = 6 K. (e) 
PL intensity 7pl after non-resonant excitation at 1.561 eV 
(T L = 4.2 K) with densities of 1 x 10 10 cm -2 (solid line) and 
3.2 x 10 10 cm -2 (dashed line). The PL intensity calculated 
from nx(t) and Tc in panels (a) and (d) is indicated for low 
and high densities, respectively, as solid and open diamonds. 



low- intensity limit, and calculations of exciton formation 
via carrier-carrier interactions in microscopic models^ 
Moreover, other nanoscale materials can be studied via 
THz probes to explore size-dependent electron-phonon 
interaction strengths and formation rates. 



B. Photoluminescence dynamics 

Next, we compare the pair kinetics with time-resolved 
PL measured in our sample, shown as lines in Fig. [7(e) for 
comparable, non-resonant excitation conditions. In stark 
contrast to the resonantly-excited case [Fig. [3(b)] , the PL 
displays a slow rise which reaches its maximum only « 
1000 ps after excitation. Such a delayed rise agrees with 



previous PL studies of exciton formation. 46 ! 50 ^^^^^ 
The fast initial exciton formation is largely absent in 
the PL, which indicates that the excitons probed by the 
THz pulses at early delays primarily populate high en- 
ergy states with momenta K ^> 0. 

There has been a long debate on whether exciton lu- 
minescence can arise in a plasma of unbound e-h pairs or 
whether it is fully explained by "bright" excitons around 
K = 0. 62 i 63 i 64 In this respect, our observations con- 
trast sharply with a recent THz study that concluded 
the absence of THz absorption (and thus of excitons) 
with simultaneous observation of exciton PL after non- 
resonant excitation^ In that study, the THz absorption 
was probed around a single wavelength assumed to co- 
incide with the ls-2p transition. Moreover, the near-IR 
exciton line was almost 10 x broader than in our sample, 
pointing to significant inhomogeneities that can result in 
carrier localization. In contrast, the present experiments 
provide full and broadband spectral information of the 
THz response, in both real and imaginary parts. Our re- 
sults, obtained in a high-quality QW sample, clearly re- 
solve the existence of a significant THz absorption peak 
at the \s-2p transition - revealing a large exciton den- 
sity - directly after non-resonant excitation, while the 
PL rises only slowly. 

Indeed, given the THz-derived densities we can test 
whether our measured PL dynamics is fully explained if 
we assume that only excitons contribute to the lumines- 
cence. Following previous work (see, e.g., Ref. I46T ). the 
intensity of the luminescent fraction can in this case be 
written as 



J PL (t) cx 



i 



-hr h /k B T c (t) 



nx(t) 
k B T c {t) 



(12) 



where Th is the near-IR homogeneous linewidth, and the 
intuitively simple approximation on the right hand side 
is valid for ffh <C ksTc- We verified that the shape of 
the PL dynamics is not affected by the approximation, 
by comparing it with the full expression in Eq. [T^with a 
density-dependent For the calculation, we use the 

experimentally-derived nx(t) [Fig. [7(a)], and the above- 
discussed Tc(t) [Fig. [7(d)] which resulted in a close de- 
scription of fx(t)- The resulting calculated PL intensity 
is shown as diamonds in Fig. [7(e). It provides a good 
representation of the luminescence rise time, underscor- 
ing the consistent agreement between the PL and THz 
signals in the above scenario. Remaining differences in 
the calculated and measured PL shapes can be explained 
in part by the lower time resolution of the PL experiment, 
and by inherent limitations of our simplified analysis that 
assumes fully thermalized exciton and free-carrier distri- 
butions. This comparison quantitatively confirms the 
sensitivity of luminescence to optically-active excitons 
around K ~ 0, leading to the predominance of relax- 
ation to low-energy states rather than exciton formation 
in the PL kinetics. 

Thus, our two-component THz analysis provides a 
fully-consistent description of exciton and free-carrier dy- 
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namics across a large set of experimental data. The ap- 
plicability of this model is underscored by: (i) the abso- 
lute density scaling, corroborated by excellent agreement 
between the absorbed photon flux and the densities ob- 
tained from model fits; with excitons and free carriers 
gauged separately, (ii) the close agreement in shape with 
measured spectra in both real and imaginary parts, {lit) 
the precise fulfillment of the intraband sum rule, (iv) 
the quantitative agreement of the THz-derived exciton 
fraction after ionization with thermodynamic equilibrium 
(Saha equation) at a consistent absolute pair density, (v) 
the likewise quantitative agreement of the THz-derived 
ionization rate with the known LO-phonon ionization fre- 
quency, (vi) the pair kinetics during exciton formation 
which follows closely a Saha equilibrium, a physically 
reasonable result assuming rapid formation at these den- 
sities, and finally as explained above (nil) the consistent 
and quantitative description of PL dynamics using the 
THz-derived pair densities. 



n=l 



N times 



FIG. 8: Multi-layer geometry consisting of a thick substrate 
(S), spacer layer (L), a thin stack of TV alternating barrier (B) 
and well (W) layers, and second spacer layer. 
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VI. CONCLUSIONS 



To conclude, we discussed optical-pump THz-probe 
studies and a detailed model analysis of the transient 
THz spectra and density kinetics of quasi-2D excitons 
and unbound e-h pairs in GaAs quantum wells. An intra- 
excitonic dielectric function was presented, whose shape 
and absolute density scaling is in excellent agreement 
with the measured THz response of resonantly-generated 
excitons. Ionization of excitons, in turn, leads to THz 
spectra that exhibit both intra-excitonic and Drude-like 
features. Here, a two-component dielectric function suc- 
cessfully describes the complex THz spectra and yields 
densities of excitons and unbound e-h pairs. Ionization is 
found to result in a quasi-equilibrium, whose exciton frac- 
tion quantitatively agrees with the Saha thermodynamic 
equilibrium - thus experimentally verifying the density 
scaling of the two-component model. 

The analysis is equally well applied to transient THz 
spectra during exciton formation, demonstrating fast ini- 
tial formation of ~ 40% excitons followed by slower pair 
binding within several 100 ps. At the longest delays, 
about 90% of the pairs are bound into excitons. The 
time-dependent exciton fraction at our densities is quan- 
titatively described by a Saha equilibrium that follows 
the e-h gas cooling dynamics. Finally, in this scenario 
a consistent agreement is found with time-resolved PL 
measured for comparison. Our study provides the basis 
for further exploration of the pair kinetics of quasi-2D e-h 
gases via transient THz spectroscopy and holds promise 
for enabling new studies of exciton physics inaccessible 
with near-IR light. 



APPENDIX A: COMPLEX THZ TRANSMISSION 
FUNCTION OF MULTIPLE QUANTUM WELLS 

In the following, we deduce the complex THz trans- 
mission function (and its pump-induced change) for a 
multiple-quantum- well structure. Figure [5] shows the ge- 
ometry, modeled as a thin multi-layer stack on a thick 
dielectric substrate. This treatment, rather than a single- 
layer approximation, becomes relevant at the lowest THz 
frequencies, since each absorbing layer adds a Fresnel 
phase shift that scales approximately inversely with fre- 
quency. We can write the complex transmission of the 
layered system using a matrix approach 66 

Ana pifn s d s 

(5n + oi 2 )n s + S 2 2 + &21 n s + 1 

where S is a matrix describing the response of the multi- 
layer stack, while ns and ds are the substrate refrac- 
tive index and thickness, respectively. This expression 
takes into account all reflections within the thin stack 
but ignores multiple reflections in the substrate, since our 
electro-optically sampled THz field trace covers only the 
first pulse replica. For the above geometry S = L • M • L, 
where the matrix 

L=( C0S(3 \ -4l si »M (A2) 
\-in L sm(3 L cos^l / 

corresponds to each of the two spacer layers and 



11 



describes the multi-layer sequence of quantum wells and 
barriers, with 



Mn = Wjv-i(a) ( cos /3b cos /?w — sin (3 B sin /3w 

n B 



M 
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-UN-2(a) 

-iUif-i(a) 



cos /3b sin /3w sin /3b cos /3w 



n w 



»B 



M 2 i = -iWjv-i(a) x 

(nB sin /3b cos /3w + nw cos /3b sin /3w) 



M22 = Un-i{o) ( cos/3 B cos/3w sin/3 B sin/5 w 

n w 



-Un-2{o) 



(A4) 



In the above, iV is the number of barrier-QW layer pairs, 
while rii is the refractive index and di is the thickness 
of layer i (with i = W, B, L denoting well, barrier, 
and spacer layer, respectively). Moreover, /3, = ^njdj 
are phase factors, and Un{x) = (1 — x 2 )^ 1 / 2 sin[(iV + 
1) arccos(x)] are the Chebyshev polynomials of the sec- 
ond kind, — with 



a = cos /3b cos /?w — r (— - + — — ) sin /3b sin f3w- (A5) 
2 nw ^b 

Approximating the formulas for our case of optically-thin 
media, -C 1, yields a « 1 and Un{o) — ► 2V + 1. 
Terms of the order of (3f and higher can be accordingly 
neglected, such that 



-i^d L n 2 h 



-i-d L 

C 

1 



(A6) 



M 



1 -iN^(d B + d w ) 

-i^(d B n| + d w n^) 1 



(A7) 

Combining Eqs. I AH IA6[ and IA71 and neglecting terms 
of the order of (^d^n-^) 2 and (^b^b) 2 or higher yields 
the complex THz transmission coefficient 



t(w) = 



4n P 



,0s 



1 + n s - e(w)A + B n s + l 



(A8) 



with 



A = i-D QW [l-i-d L (l + n s )} 
c c 

fl = -i-(n s Aot + 2d L n 2 + Nd B n B ) 
c 

where Dtot = 2g?l + N(d B + dw) is the total multi-layer 
thickness and Dqw = A^dw is the aggregate thickness 
of quantum well material. In the above, e(ui) = n^y(w) 
is the complex-valued quantum-well dielectric function. 
Note that in the limit c?l, d B — > 0, the above expression 
reduces to the transmission of a single optically-thin, ab- 
sorbing layer on a dielectric substrate. 



In equilibrium, Eq. IA8I directly connects the static 
dielectric function e(co) to the experimentally-accessible 
complex transmission coefficient given by t(u>) — 
E(uS) / E- m {u)) . Here, E- ln {u>) and E(u>) are the incom- 
ing and transmitted fields, respectively. Likewise, in 
the photoexcited state the modified dielectric function 
e(uj) + Ae(uj) is linked to t*(uj) = [E(tv) + AE(uj)}/E in (u}), 
where AE(u>) is the pump-induced field change. The ra- 
tio of these transmission coefficients is then given by 



*M = 

t*(u) E{uj) + AE(uj)' 



(A9) 



Combined with Eq. IA8[ this yields an analytical expres- 
sion for the pump-induced change in the quantum well 
THz dielectric function, 



Ae(w) 



E{uj) + AE(cj) 



l + n s + B 
A 



(A10) 



The above multi-layer expression accounts for the effects 
of Fresnel phase shifts at each interface. For our specific 
structure, we have N = 10, d^ — 500 nm, d B = 10 nm, 
dw = 14 nm, ng =3.1, = n B =3.4, and (in equilib- 
rium) nw =3.6. While Eqs. lA8l and lA10l are sufficient ap- 
proximations for most conditions, the dielectric function 
change Ae(<x>) in the high-density regime can be obtained 
via numerical solution of Eqs. lAll and IA91 



APPENDIX B: INTRA-EXCITON TERAHERTZ 
DIELECTRIC RESPONSE 

Below, we provide a detailed derivation of the THz di- 
electric response in Eq. [6] due to transitions from the 
Is exciton state into higher bound states and into the 
continuum. Particular attention is paid to scaling the 
response in absolute units. The bound 2D exciton nor- 
malized wave functions ar o 36 ' 67 i 68 



S 2 (n-\m\)\ 
a 2 (n+±) 3 ' ((n+|m|)!) 3 



xp 1 



I p/2 L 2\m\ ,v 
n+|m| / 



7^ 



(Bl) 



the main quantum 

\m\ < n, LP{p) = 

v)\(p + v)\\ are the as- 



where n = 0,1,2... indicates 
number, m is an integer with 

sociated Laguerre polynomials, and p = 2r/((n+ l/2)a). 
The Bohr radius a and binding energies E n scale as in 
Eqs. 2] and [5] Following Ref. [3y, the finite well size is 
taken into account by rescaling the Coulomb potential 
with a parameter A, in order to recover a realistic bind- 
ing energy. The quantum number n is enumerated start- 
ing from zero, while we will colloquially call the ground 
state "Is" (corresponding to n = 0,m = 0), and the 
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higher bound states accordingly "n + 1" levels, i.e., the 
"2p" level corresponds ton=l,m = ±l. 

The wave functions of unbound e-h pairs in the con- 
tinuum are, in turn, 



i>k, m {r, 4>) 



(2fcr)l r 



(2 H )!" F l H + ^ + i ;2! " " 1;2//rr 



xe 



2A- 



l + e 



I 7 "-! 

n 

i=i L 



C? -2 )2 + ^ 



(B2) 



where F(a;b;z) denotes the confluent hypergeometric 
function, and the product is replaced by unity for m — 0. 
These functions are normalized per unit momentum, 
i^k' m'\' l l ; k,m) = S(k — k')5{m — rn'). The scalar wave 
number k is defined via Uk = ^/2pE(k), where E(k) is 
the unbound e-h pair kinetic energy relative to the band 
edge. 

Using Fermi's golden rule, the contribution to the di- 
electric function e(u>) for excitons in the Is ground state 
ipo.o is given by 



Ae(w) 



/i 



([ 



g +g(fc) l2 



iu;T bh 

- dk) , (B3) 



which represents the sum over all final states ipf with 
oscillator strengths 



/ls,f 



2// 



{E { -Eo) \(^\x\^ 0fi }\' 



(B4) 



In the above, nx is the Is exciton sheet density, \i is the 
reduced e-h pair effective mass, and Tbb, Tb c are phe- 
nomenological broadening parameters for bound-bound 
and bound-continuum transitions, respectively. 

Consider first transitions to p-like bound exciton 
states. For an in-plane THz field linearly polarized along 
x = r cos(^), the dipole matrix elements for each of two- 
fold degenerate final states (to = ±1) are 



(lpn,±l\x\lpO,o) = 



a I (n+|)3(n- l)l ;< 



((n + 1)!) 3 



(B5) 

,3 e -(„ +1 )„ L 2 (p) dp 



The latter integral can be written ana- 
lyticaUy as p 3 e -^PL 2 n+1 (p) dp 



(n - l)!n n (2n + l)(n + l)"" -1 , 
tor strengths 



yielding the oscilla- 



(B6) 



when including both to = ±1 final states. The matrix 
elements for transitions from Is into p-like continuum 
states are, in turn, given by 



2k 2k , I 1 , 

< Vfc,±iR#o,o >= — \ — (t + ~ 2T2 ) 

a V l + er~ 4 a 2 fc 2 

-«"> 3 F[- + — ;3;2<Ar](/r 
V 2 ak 



x e 
Jo 



(B7) 



The above integral can be written as 
J~ e -2r/o-ifcr r 3 F (| + A. . 3. 2ifcr)dr = - 8to 4 (1 + 

2ak/(2i — ak))~^~^ (— 2« + afc) -5 , yielding the oscillator 
strengths per unit momentum 



2ua 2 256a 2 fc(2i±4)"^ 

/ ls ( fc ) = (£ + E(k)) 2J 2i ~ afc 

V ' h 2 V W7 (l + e-^)(4 + a 2 P) 4 



(B8) 



as the sum of transitions to to = ±1. With the above, 
we obtain the dielectric response Ae(w) in Eq. [6] arising 
from intra-excitonic transitions between the Is state and 
all higher bound states and the continuum. 

Finally, we can verify the expected smooth transition 
of the dielectric function at the edge between bound and 
continuum transitions, in analogy to the intra-atomic ab- 
sorption of hydrogen^ The theory used to derive the 2D 
exciton wave functions links the quantum number n and 
wave number k via = 1/ {ia^k). For bound states, n 
is integer and real, while for continuum states it is a con- 
tinuous and imaginary number. To compare the two com- 
ponents in Eq.[51 the oscillator strengths must be normal- 
ized to a constant "quantum number" interval An = 1, 
i.e., the fc-dependent dipole moment must be multi- 
plied by (dk/dn)An = —l/(iao(n + \) 2 )- Substituting 
l/(wofc) — ► n +5 then transforms the bound-continuum 
expression into one that differs from the bound-bound ex- 
pression only by a factor of (1 + e - 27r l™+2 1) -1 . Since the 
latter approaches unity for n — > 00, both expressions ana- 
lytically match up smoothly at the continuum edge. This 
fact is illustrated in Fig. [^g) by numerical simulation. 
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